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We study the collapse of a free scalar field in the Brans- 
Dicke model of gravity. At the critical point of black hole 
formation, the model admits two distinctive solutions depen- 
dent on the value of the coupling parameter. We find one 
solution to be discretely self-similar and the other to exhibit 
continuous self-similarity. 



Studies of black hole formation from the gravitational 
collapse of a massless scalar field have revealed inter- 
esting nonlinear phenomena at the threshold of black 
hole formation . These studies have shown that Ein- 
stein's field equations possess solutions which occur pre- 
cisely at the black hole threshold and which are univer- 
sal with respect to the initial conditions of the evolution. 
More specifically, for any type of initial field configura- 
tion whose energy is parameterized by some parameter, 
p, the critical solution occurs at a value oi p = p* such 
that for all p < p* no black hole is formed, and for all 
p > p* a black hole is necessarily formed. This critical 
solution, whether obtained with an initial pulse shape 
such as tanh or a Gaussian pulse, is identical, erasing all 
detail of the initial field configuration. 

Though universal with respect to initial conditions, the 
critical solution is dependent on the specific matter model 
involved. In the case of a real scalar field 0| , a discretely 
self-similar solution (DSS) is found, characterized by an 
echoing exponent A. In other words, were an observer 
to take a snap-shot of the solution at some time t, he 
would find the same picture as when he zoomed in to a 
spatial scale exp(A) smaller than the original at a time 
t + exp(A) later. 

In contrast to this DSS solution, other researchers, 
working in an axion/dilaton model, have found that the 
equations possess a continuously self-similar (CSS) so- 
lution H]. Because they found this solution by assum- 
ing continuous self-similarity and solving the appropri- 
ate ordinary differential equations, they could not show 
whether this CSS solution is indeed a critical solution. 

We find that a free real scalar field coupled to Brans- 
Dicke gravity contains two distinct dynamic critical so- 
lutions. As a special case, the model includes the real 
scalar field in general relativity and recovers the DSS so- 
lution as in Further, this model is sufficiently general 
that it contains the model studied in j^] as another spe- 
cial case. For this case, we find that the CSS solution 
is an attracting critical solution. Hence we present the 



novel result that for a single matter model, adjustment 
of a coupling parameter transitions between two unique, 
dynamic, attracting critical solutions. Because these two 
solutions are both dynamic, the model is quite different 
from the Yang Mills model studied in Q . 

Subsequent to our study, Hirschmann and Eardley, 
working in an even more general model, the non-linear 
sigma model, which includes ours, carry-out a pertur- 
bation analysis and confirm a change in stability near 
the value we find for the transition coupling parameter 
Further, from the eigenvalues of the unstable modes, 
they have been able to compute mass-scaling exponents. 
Their results concur with those we find from our numer- 
ical evolutions. 

We work in spherical symmetry with the metric 



ds^ = -a(r, tfdt^ + a{r, tfdr^ + r^d^^ 



(1) 



where a{r,t) represents the lapse function in the 3+1 
formalism and r measures proper surface area. 

The Brans-Dicke model is described by the field equa- 
tions 
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where l/0(r, t) represents the freedom of the conventional 
gravitational constant to vary ||^ . The total stress-energy 
tensor consists of two terms 
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where T^^ represents the energy associated with the 
Brans-Dicke field (j> and is the conventional ten- 

sor associated with matter sources Q . For this study our 
sole matter source is a free massless minimally coupled 
scalar field ipir, t) governed by 
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and whose associated stress-energy is 
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The Brans-Dicke field satisfies the generally covariant 
wave equation 



0(h = AttXT" 
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where A, a constant, represents the strength of the cou- 
pling between the Brans-Dicke field and matter Its 
associated stress-energy tensor is 
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where 
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The equations described above are said to be in the 
Brans-Dicke frame where masses are constant but iner- 
tial forces depend on the distribution of mass in the uni- 
verse. However, it is possible to transform to a conformal 
frame in which the geometry is described by Einstein's 
field equations with vanishing second derivatives of 4>. In 
this frame, the Einstein frame, masses vary with time, 
but the gravitational constant is indeed constant. 
We achieve this conformal transformation via 
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after which we have the equations (now expressed in the 
Einstein frame) 
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We define auxiliary variables in terms of the derivatives 
of the scalar fields 
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so that the wave equations result in 
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The only other necessary conditions come from the field 
equations, which, in the Einstein frame, are simply Ein- 
stein's field equations, G^j/ = SttT^;^. In accordance with 
the 3+1 formalism, we have the Hamiltonian constraint 
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and the polar slicing condition 
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which enable us to solve for the geometry in terms of the 
two sources, ip and ^. These equations suffice to evolve 
both the fields ip{r, t) and ^(r, t), and the geometric vari- 
ables a{r,t) and a{r,t) |^]. 

To show that the model found in Q is a special case 
of our model, we compare our Lagrangian 
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with that of 
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defined in terms of the axion, a, and the dilaton, (p. Com- 
paring Eqs. ( p9| ) and (pO|), we see a correspondence be- 
tween the two models with a trivial rescaling of the fields 

^ = -40 V = ^^a A = 8. (21) 



We have found the critical solutions for a variety of ini- 
tial data. Specifically, we input the initial configuration 
of the two fields, and specify the value of A. The space 
of initial configurations is schematically represented in 
Fig. 0. 
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FIG. 1. Schematic of the configuration space. The hori- 
zontal axis represents the mixing ratio between the free scalar 
field and the Brans-Dicke field. The darkened triangles rep- 
resent the locations of the profiles displayed in Fig. 2. The 
darkened circle represents the location of the real scalar field 
in general relativity studied in [1]. For A < 2/3, only the DSS 
solution is the attractor. Above A ~ 2/3 the CSS solution 
attracts whenever both fields are initially present. 
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We observe for an initially vanishing scalar field that 
Eqs. 1^ |l^ describe the real scalar field case studied 
in Consistent with this observation, our results re- 
cover the same DSS solution found for the real scalar 
field case. The equivalence between this model with 
■!/)(r, t) = and that of the real scalar field is shown in 
Fig. as the vertical line extending through the middle 
of the graph. 

When A 0+, Weinberg shows that the Brans-Dicke 
model goes over to general relativity. Hence, for the gen- 
eral situation in which both fields are present {tp/^ 0), 
we expect to recover the results from general relativity. 
We do recover the general relativity result, that being the 
DSS solution. As shown in Fig. 0, the critical solution 
is discrete for generic initial data as A is increased up to 
Aw 2/3. 

Shown in Fig. || for A = 0.09, we have verified that 
this is the same DSS solution obtained for the real scalar 
field in general relativity In Table |, we show the 
computed values of the echoing exponent A. These val- 
ues correspond to that found in . 




FIG. 2. Demonstration of regime in which the solution 
transitions from the DSS to the CSS. We show here the 
mid-section of the seven solutions for various A. For A = 0.09 
the solution is clearly the DSS, however the next solution 
demonstrates that the echoes get damped as one moves to- 
wards the origin. Eventually the solutions become the CSS. 

Around A « 2/3, a remarkable transition occurs in the 
critical solution. As one increases A in this region, the 
echos displayed by the critical solution are damped by a 
decreasing envelope as shown in Fig. |^. 

At A = 8, we recover, as expected from Eq. (pi]), the 
CSS solution found in 0. In Fig. ^we demonstrate that 
the solution found by by demanding continuous self- 
similarity is indeed the attracting critical solution. Here 
we show that by a trivial rescaling of the fields at one 
time slice, our solution is identical to theirs. 
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FIG. 3. Demonstration of equivalence between the CSS so- 
lution obtained from [2] and our critical solution obtained 
with A = 8. The solid line shows the metric function a(r, t) 
versus Inr provided by Eardley. The crosses denote data 
points from our solution. Four time profiles are shown with 
T = ln(T — T*), where T is the central proper time of the 
slice and T* is the critical time of collapse. The Eardley 
and Hirschmann solution is scaled to match our profile at 
T — —3.3. The congruence at other times displays the equiv- 
alence of the two solutions. 
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FIG. 4. Black hole mass scaling exponents for various A. 
The dotted line displays the values obtained for the non-linear 
sigma model in [4]. The open pentagons represent the scaling 
exponent obtained by least-squares fits using our numerical 
results. The errorbars represent a range of three standard 
deviations. 
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FIG. 5. Illustration of the power-law mass scaling relation. 
The markers display the mass obtained for the normalized 
distance from criticality. The lines designate the least-squares 
fit line with slope 7. 

Both the critical solutions exhibit mass-scaling in the 
supercritical regime. Specifically, for some family of so- 
lutions where p* represents the critical value of a param- 
eter, the masses of the black holes formed in the regime 
where p > p* follow 
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where 7 depends on A and c is a family dependent con- 
stant. Fig. [s] shows four power law fits and the associated 
7's. 

In keeping with the correspondence between this model 
for very small A and that of general relativity, we find 
7 — 0.37, matching that found in Likewise, we find 
agreement between our values of 7 and those found by 
perturbation analysis in We display both these sets 
of values in Fig. ^. 

The appearance of these two disparate solutions leads 
one to examine the transition in A-space from the DSS 
to the CSS. Bracketing solutions have shown that around 
A = 2/3 the transition occurs (see Fig. As A is in- 
creased around this transition value, an envelope damp- 
ens the discrete echos into the smoothly continuous self- 
similar solution. Perturbation results in the non-linear 
sigma model confirm a change in stability of the CSS 
solution near A = 2/3 

Further parameter surveys are needed to further spec- 
ify the transition point between the two self-similar solu- 
tions. We also anticipate interesting negative mass solu- 
tions for A < 0. However, our studies have clearly shown 
the richness of the solution space for even a simple, two- 
scalar field, one-dimcnsional problem such as this one. 
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TABLE I. Mass scaling exponents 7 and the spatial scaling 
exponents A for the various discretely self-similar solutions 
found. 
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0.001 


1000 


0.375 


0.01 


100 


0.379 


0.087 


10 


0.374 


0.4 


1.0 


0.374 


0.5 


0.5 


0.378 


0.57 


0.25 


0.373 


0.615 


0.125 


0.373 


0.67 
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0.375 


0.73 


-0.125 


0.372 


0.8 


-0.25 


0.371 


1.0 


-0.50 


0.364 


1.33 


-0.75 


0.362 


2.0 


-1.00 


0.348 


4.0 


-1.25 


0.314 


8.0 


-1.375 


0.263 


10 


-1.40 


0.245 


20 


-1.45 


0.205 



A 



3.447 
3.447 
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